Real Space Hartree-Fock Configuration Interaction Method For Complex Lateral 

Quantum Dot Molecules 
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We present unrestricted Hartree Fock method coupled with configuration interaction (CI) method 
(URHF-CI) suitable for the calculation of ground and excited states of large number of electrons 
localized by complex gate potentials in quasi-two-dimensional quantum dot molecules. The method 
employs real space finite difference method, incorporating strong magnetic field, for the calculating 
of single particle states. The Hartree-Fock method is employed for the calculation of direct and 
exchange interaction contribution to the ground state energy. The effects of correlations are included 
in energies and directly in the many-particle wavefunctions via configuration interaction (CI) method 
using a limited set of excitations above the Fermi level. The URHF-CI method and its performance 
are illustrated on the example of ten electrons confined in a two-dimensional quantum dot molecule. 

PACS numbers: 73.43.Lp 



I. INTRODUCTION 

There is currently interest in localizing electrons in 
quantum dots and utilizing their spin for quantum 
information processingiiSi^. Alternatively, we can view 
such devices as highly tunable quantum dot molecules, 
with large numbers of electrons or qubits, with gates 
and magnetic field used to perform quantum operations. 
Even with a relatively small number of bits, to under- 
stand working of such devices one needs a computational 
tool suitable for the calculation of wavefunctions of 
at least tens of electrons confined by complex gate 
potentials and in magnetic field. The wavefunctions 
must capture entanglement, or in condensed matter 
language, correlations. The configuration interaction 
(CI) method is a suitable candidate modified here for the 
problem of quantum dot molecules. In the configuration 
interaction (CI) method the Hamiltonian of an interact- 
ing system is calculated in the basis of a finite number of 
many-electron configurations and diagonalized exactly. 
Exact diagonalization (ED) has been an important 
tool in the context of quantum chemistry of atoms and 
molecules*. In condensed matter physics ED has been 
used to investigate the electronic and optical proper- 
ties of artificial atoms and molecules in quantum dots 
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and finite size quantum Hall systemsi^SiSSi^SiiS The di- 
rect access to electron correlations in the many-particle 
wave functions calculated by CI makes CI useful for 
the purpose of quantitative description of the electron 
entanglement in the context of solid state quantum 
information architectures. The drawback of the method 
is its slow convergence and applicability to a small 
number of electrons. While mean-field theories such as 
Hartree-Fock (HF )r-b . 42 ^ ^^^-^^ density-functional theory 
(DFT)43,44i45,46,47,48,49 applied to large elec- 

tron numbers, they neglect electron correlations, or 
entanglement, in the wavefunction. 

To include electron correlations explicitly in the many- 
particle wavefunctions one may follow quantum chem- 



istry methods which start with HF or DFT calculations of 
effective single particle orbitals, and use them in the con- 
struction of a finite number of many-electron configura- 
tions as input in configuration interaction calculations^. 
These ideas have been already explored for the calcula- 
tion of electronic properties of quantum dots with large 
electron numbers^^ and quantum dot molecules 1^2*2^ 

The objective of this work is to present an efficient HF- 
CI method: (a) based on real space calculations of single 
particle states and hence suitable for complex quantum 
dot potentials, (b) combining complex potentials with 
strong magnetic field, and (c) capable of treating corre- 
lations (entanglement) in electron wavefunction for large 
number of electrons. The URHF-CI method and its per- 
formance are illustrated on the example of ten electron 
confined in a two-dimensional quantum dot molecule. 

The paper is organized as follows: In section ^1 we 
present the Hamiltonian of electrons confined in the lat- 
eral gated quantum dot molecule. The calculation of the 
single particle spectrum, including the details of the im- 
plementation of the localized gauge, is presented in sec- 
tion mil In section IIVI our many electron approach, in- 
cluding the unrestricted Hartree-Fock approximation, is 
introduced. The central results of this study, the cal- 
culation of the ground state energy and excitation gap 
.29.iiixa]c8aiferSiii:M.i3fr i interaction method using unrestricted 



Hartree-Fock URHF-CI basis and non-interacting single 
particle SP-CI basis, are presented in section The 
paper is summarized in section IVTl 



II. HAMILTONIAN 

We consider N electrons with effective mass to*, mov- 
ing on a plane {x,y) defined by the GaAs/GaAlAs in- 
terface, with complicated confining potential V{r) as a 
function of electron position r, created by metallic gates 
on the surface of the hcterostructure. The magnetic field, 
defined by vector potential A, is perpendicular to the 2D 
plane. The effective TV-electron quantum dot molecule 
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FIG. 1: discrete points in left (L) and right (R) dot regions. 



Hamiltonian in the presence of an external magnetic field 
can be written as: 



N 



Zm* \ t c 



(1) 



where e is electron charge, c is velocity of light, and e is 
the host semiconductor dielectric constant. The Zeeman 
splitting (very small for GaAs) is neglected here. The 
first term in Eq. is the sum of single particle Hamil- 
tonians and the second term is the Coulomb repulsion. 
The electron-electron repulsion is screened by the dielec- 
tric constant, which for GaAs is e « 12.4. 

In what follows, we present the results in effective 
atomic units with length unit in effective Bohr radii 
ag = (?i^/2m*)/(e^/2e), and energy unit in effective Ry- 
dberg Ry* ~ e^/(2eag). Using material parameters for 
GaAs with m* k, O.OGTtoq (mo is the bare electron mass) 
Ry* = 5.93meV, and = 9.79nm. 



III. SINGLE PARTICLE SPECTRUM IN A 
MAGNETIC FIELD 

The single particle Schrodinger equation is given by 
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ifi{r) = £(^(f), (2) 



where A = {B /2){—y, x) is the vector potential seen 
by electrons, V{r) is the confining potential, and B is 
the magnetic field. For a single quantum dot the po- 
tential is often well approximated by an analytical form 
V{f) — l/2mcj^r^. For quantum dot molecules the po- 
tential will have several minima in which to localize elec- 
trons, and the most convenient way to solve the single 
particle problem is by discretizing the Schrodinger equa- 
tion. As an example of a complicated molecular quan- 
tum dot potential we study a double Gaussian potential 



R exp[ ^ 

with Vl and Vr measuring the 
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depth of the left and right dot, and Vp controlling the 
height of the tunneling barrier. The centers of the two 
dots are separated by 2a. To describe the typical cou- 
pled quantum dots we use the parameters Vl = Vr = 
-10, a = 2, A = 2.5, and Ap^ = 0.3, Ap.y = 2.5, in 
effective atomic units. Vp which controls the central po- 
tential barrier mimics the plunger gate strength and is 
varied between zero and 10i?y*, independent of the loca- 
tions of the quantum dots. 

The choice of gauge A plays significant role in improv- 
ing the numerical accuracy of single particle spectrum. 
We use a gauge field which is adopted to the geometry of 
the confining potential. For a double dot molecule with 
two minima, we divide the plane of coupled quantum dots 
into left and right domains, as shown in Fig. ^ In or- 
der to end up with a well defined left and right quantum 
dots with increasing barrier height, i.e. with dots with 
gauge centered in the origin, we define left (right) wave- 
functions and Hamiltonians, and carry out gauge trans- 
formation on both Hamiltonians and wavefunctions. The 
resulting vector potentials localized in centers of each dot 
are Al — {B/2){~y, x + a) for the left dot { x <0), and 
Ar = {B/2){-y,x - a) for the right dot {x > 0). Intro- 
ducing corresponding wavefunctions in the left (right) dot 
'P^.j = '^f^j exp(-|-iWca2//4) (tp^j = ip^jexp{-iuJcay/4.)) 
allows us to write an explicit form of the discrete form 
of Schrodinger equation. In the left-dot the Schrodinger 
equation can be obtained by expanding (f^ around the 
point (z,j) 



(3) 



where Xi — ih, and yj — jh, and h is the grid spacing. A 
similar equation holds for the wavefunction in the right 
dot. 

As it is illustrated in Fig. adjacent to the inter-dot 
boundary, the point (i,j) in L, is connected to the point 
[i + l,j) in R, via Schrodinger equation. In this case 
'f'i+i^j in Eq.© is already in the right dot and is not 
known and must be replaced by 'fi+i,j exp(— ja;cay/2). 

One should note that this separation of gauge fields 
into left and right leads to a uniform magnetic field B — 
V X A in both half-planes, but it produces an artifical 
discontinuity in B along a; = 0. To avoid this unphysical 
discontinuity, in finite difference method, the grids along 
the a;-axis have been set slightly away from x = such 
that X = 0-line has been excluded from the real space 
finite difference calculation. 

The single particle spectrum and (pj is calculated 
accurately on a finite mesh of x Ny (102 x 51 mesh 
points are used in what follows) using conjugated gra- 
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FIG. 2: Double dot single particle spectrum vs. cyclotron 
energy. Several crossing between Landau levels is visible at 
TiLOc = lARy*. 



in a self-consistent HF field. Similar to single parti- 
cle configuration-interaction (SP-CI) method, the HF or- 
bitals can be used as a basis of the interacting Hamilto- 
nian in second quantization, and subsequently in the CI 
calculation. In this method electron-electron interactions 
are included in two steps: direct and exchange interac- 
tion using Hartree-Fock approximation, and correlations 
using HF basis in the configuration interaction method. 



A. Unrestricted Hartree-Fock Approximation 
(URHFA) 

Hartree-Fock approximation is a mean field approach 
to many body systems which accounts for the direct and 
exchange Coulomb interactions. Combining HFA with 
more sophisticated many body methods (such as CI) al- 
lows to isolate the effect introduced by correlations. The 
Hartree-Fock ground state (GS) of the electrons with 
given Sz = {N^ — iV|)/2 is a single Slater determinant 



dient methods. The magnetic field dependence of single 
particle spectrum calculated by this method, is shown in 
Fig. 121 At zero magnetic field Fig. |2shows the formation 
of hybridized S, P, and D shells. In high magnetic field 
we observe the formation of shells of closely spaced pairs 
of levels (with opposite parity). In intermediate magnetic 
field states with opposite parity cross and states with the 
same parity anticross. The crossing of many levels at e.g. 
hojc = lARy*, marked with dashed line in Fig|21 leads 
to the most computationally challenging many-electron 
problem. 



IV. MANY ELECTRON SPECTRUM 

The single-particle (SP) states calculated in previous 
section can be used as a basis in configuration-interaction 
(CI) calculation. Denoting the creation (annihilation) 
operators for electron in non- interacting SP state \aa) 
by cjj^ {caa): the Hamiltonian of an interacting system 
in second quantization can be written as 



where the first term is the single par- 
ticle Hamiltonian, and VaaSa' ,va = 
j dr ^ d?(pl^{r)(p}^^,{?)^^^^ip^,^,{?)(p^a{f), is the 

two-body Coulomb matrix element. 

Alternatively we may use the single particle states to 
first construct the Hartree-Fock orbitals. In this scheme 
electrons are treated as independent particles moving 



l*> = "^/^ E("^)^'^i'^i {rpi)\ai) . . . ipNaM {rpN)\<yN), 

(5) 

with P the permutation operator. The HF orbitals \ fia) 
which describe the state of a dressed quasi-particle in the 
quantum dot molecule can be determined by minimizing 
the HF energy Ehf = with respect to ^Piair) 

a' 3 = 1 J 

N„, 



a' 1 = 1 •' 



(6) 



~ 2jtF(7^ f^) non-interacting single par- 

ticle Hamiltonian, and V{r) — e^/er is the Coulomb 
interaction. This equation is a set of coupled equa- 
tions, for spin up and spin down states. To find numer- 
ical solutions of the Hartree-Fock equations we expand 
the HF orbitals \fia) in terms of single particle states 

|<^a): W'lccr) ^ Y.aLiO'al\'f>a,cr)- This transforms the 
HF equation, Eq.®, to the self-consistent Pople-Nesbet 
equations^: 

JVj Ni N^, 

E{^~Am+ E ^Ma/57EE«:l''«S' 

7=1 a,/3=l ct' j = l 

- w E E = (7) 

a' j = l 

where Vap^v are Coulomb matrix elements calculated us- 
ing non-interacting single particle states. 




FIG. 3: URHF eigen-energies vs. cyclotron energy for 10 
electrons with Sz = 0. The HF energy gap between HOMO 
(highest occupied molecular orbital), and LUMO (lowest un- 
occupied molecular orbital) is clearly visible. A comparison 
between SP-CI and URHF-CI calculation will be presented at 
hiJc = lARy* , and it is shown by the dash line. 



This procedure results in Ni HF states. The TV-lowest 
energy states form a Slater determinant occupied by HF 
(quasi) electrons corresponding to HF ground state. The 
rest of orbitals with higher energies are outside of the HF 
Slater determinant (unoccupied states). 

The calculated HF eigen-energies for the = 10 elec- 
trons with 5"^ = as a function of magnetic field are 
shown in Fig. |2| Comparing HF spectrum (Fig. |2Jl with 
the single particle spectrum (Fig. one observes that a 
HF gap developed at the Fermi level, between the high- 
est occupied molecular state, and the lowest unoccupied 
molecular state, and the Landau level crossing between 
single particle Landau levels has been shifted to lower 
magnetic fields. 

Fig. 21 illustrates the convergence of HF energy with 
respect to the number Ni of single particle orbitals used 
in HF expansion. 

In the following we study the effect of correlations by 
including quasielectron-quasihole excitations in the HF 
wave functions using CI method. 



V. CONFIGURATION-INTERACTION 
METHOD 

In the configuration interaction method the Hamilto- 
nian of an interacting system is calculated in the basis 
of finite number of many-electron configurations. The 
total number of configurations (or Slater determinants 
participating in CI calculation) is determined by Nc = 
[Ns\/N^\iN, - N^y.][N,\/Ni\{Ns - iV|)!]. Here 7V| and 
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FIG. 5: URHF electron-hole excitation energies Ec = 
{*Afcl-ff|*A/c) for = 13 corresponding to Nc = 1656369 
is plotted vs. Mc- Here 1 < Mc < Nc is an integer number. 
I^Pmc) is the Slater determinant of the configuration Mcth, 
formed by HF basis. Ec = Ehf (the HF ground state) if 
Mc = 1. The URHF density of states (dMc/dEc) is shown 
in inset. Its peak is below the HF energy mean value. 



Ni are the number of spin up and spin down electrons. 
This Hamiltonian is either diagonalized exactly for small 
systems or low energy eigenvalues and eigenstates are 
extracted approximately for very large number of config- 
urations. In URHF-CI (SP-CI), URHF (single parti- 
cle) orbitals are used for constructing the Slater determi- 
nants. By removing electrons from the occupied URHF 
state obtained by minimizing the total Hartree-Fock en- 
ergy and putting them onto an unoccupied URHF state, 
one can construct a number of configurations correspond- 
ing to electron-hole excitations. These excitations con- 
tribute to the many body wave functions as correlations. 
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FIG. 6: The URHF excitation energies Ec corresponding to 
the Mcth diagonal matrix elements in H (see EallO^. are 
shown for various configurations. Due to the finite band width 
W of the HP spectrum, they approximately follow a universal 
functional form. 



Denoting the creation (annihilation) operators for 
URHF quasi-particles by cj (q) with the index i repre- 
senting the combined spin-orbit quantum numbers, the 
many body Hamiltonian of the interacting system in the 
URHF basis can be written as: 



'ijklclc^CkCh 



(8) 



ijkl 



where Vijki are the Coulomb matrix elements in the 
URHF basis. Our method of computing Coulomb ma- 
trix elements is presented in the appendix. Here 



{i\T\j)=eA,-{t\VH + Vx\3), 



(9) 



where are the URHF eigenenegies, and Vh and Vx are 
the Hartree and exchange operators. The Hamiltonian 
matrix is constructed in the basis of configurations with 
definite Sz, and diagonalized using conjugated gradient 
methods. As the size of URHF basis goes to infinity, the 
method becomes exact. In practice, one may, however, 
be able to obtain accurate results using finite number of 
basis in CI calculation. 

To introduce a systematic method for constructing the 
Hamiltonian, we select a class of configurations which 
have the highest contribution to the ground state wave- 
function. Because the main contribution to the ground 
state energy of many electron system comes from the 
Slater determinants formed by the lowest energy HF or- 
bitals, we accept Slater determinants whose HF energies 
are below an energy cut-off Ec- The number of such 
Slater determinants is finite, and is given by Mq- 

The many body Hamiltonian matrix constructed in 
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FIG. 7: The ground state energy and the first excited state 
energy in URHF-CI method is shown. 



this way is the following: 



H21 



H = 



H12 
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Q 



(10) 



Here Hu < H22 < ■ ■ ■ < HmcMc < ■ ■ ■ < H^^n^ 
are the HF matrix elements with Hn = Ehf, and 
HmcMc — (the Coulomb interactions in the diag- 
onal elements account for the direct and exchange inter- 
actions). 
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FIG. 8: The dependence of URHF-CI and SP-CI ground 
state energies on Ns. The pre-calculated direct and exchange 
energies in URHF-CI improves the ground state energy and 
the convergence rate. Up to Ns = 10 the SP-CI and URHF-CI 
Hamiltonian have been diagonalized exactly. Above Ns — 10 
the variational ground state energies are calculated by diago- 
nalizing Heff with Mc = 63504. 



To test this method numerically we set the cyclotron 
energy to hcuc — lARy*, which corresponds to the cross- 
ing and mixing of many orbitals, as shown in Fig. |21 
Fig. [S] shows the evolution of URHF electron-hole exci- 
tation spectrum as a function of cut-off Mc- The energy 
bandwidth W = H^^j^^ — Ehf is finite. As it is illus- 
trated in the inset of Fig. [3 the peak of URHF density of 
states (DOS) is in the middle of the band. As it is shown 
in Fig. El the diagonal elements of the Hamiltonian H 
(URHF eigen-energies) can be represented approximately 
by Ec — Ehf + Wf{Mc/Nc) where f{x) is a universal 
function. By increasing Nc, the total number of URHF 
eigen-states Mc (corresponding to given Ec) increases. 
In this case, it is possible to find a lower energy varia- 
tional wave function if the energy cutoff of H^s is set to 
a given Ec- 

We now turn to investigate the convergence rate of the 
calculated ground state energy as a function of cutoff 
Nc and Mc- In Fig. □ the energies, measured from 
the HF energy, of the ground and the first excited state 
for = 10 electrons obtained using URHF basis with 
Ns — 10 and with Ng = 15 orbitals. We note that for 
A^ = 10 and A^ = 10 the number of configurations Nc = 
63504 while for N, = 15 and N = 10, Nc = 9018009. 
Hence for Ng — 10 the Hamiltonian can be diagonalized 
exactly and the ground and excited states are known up 
to Mc = Nc- This is not the case for Ng = 15 where we 
were able to extract the ground and excited states up to 
Mc = O.OOGAc, and exact energies are not known. 

However, in both cases the energies fall off rapidly and 
very quickly saturate as a function of Mc/Nc- Hence 
it appears sufhcient to use only a fraction of low energy 
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FIG. 9: The exchange energy, the spin singlet-triplet energy 
gap in URHF-CI method is shown. At huc = l-4Ry* the spin 
singlet state is the ground state. 



configurations to construct the effective Hamiltonian Hcs 
in order to achieve satisfactory convergence. 

Given the convergence criteria established above, we 
now discuss the advantages of URHF-CI versus much 
easier to use SP-CI, a central result of this work. Fig. 
|S1 shows the dependence of URHF-CI and SP-CI ground 
state energy on the number of single particle orbitals Ng 
for A = 10 electron droplet at magnetic field correspond- 
ing to huic = 1.4i?y*. Up to Ns = 10 the SP-CI and 
URHF-CI Hamiltonian have been diagonalized exactly. 
Above Ng = 10 the variational ground state energies are 
calculated by diagonalizing Hcs with Mc ~ 63504. The 
lowest number of single particle orbitals populated by 
A = 10 electrons with S", = is A. = 5. For A, = 5 in 
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FIG. 10: A comparison between the exchange coupling J 
calculated in SP-CI and URHF-CI is shown. Up to Ns = 10 
the SP-CI and URHF-CI Hamiltonian have been diagonalized 
exactly. Above TVs = 10 a variational calculation has been 
used by diagonalizing HcH with Ale = 63504. 



SP-CI electrons populate the lowest five single particle or- 
bitals while in URHF electrons populate the HF orbitals, 
which minimize not only single particle energy but also 
direct and exchange energy. Hence the starting energy of 
initial single configuration in URHF-CI has significantly 
lower energy compared to SP-CI. As the number of avail- 
able states Ns, and hence the number of electron- hole 
excitations (configurations) increases, the ground state 
energy obtained in URHF-CI decreases very slowly. It 
starts with Hartree-Fock value of E ~ ~29.13Ry* for 
Ns = 5 and ends up with E = -29.6Ry* for N^ = 15. 
This gives our best estimate of total correlation energy 
of Ec — 0A7Ry*. We find the correlation energy to be 
only two percent of Coulomb energy. 

By contrast with URHF-CI the SP-CI calculations us- 
ing the single particle basis converge very slowly. The 
slow convergence can be understood in terms of large di- 
rect and exchange energy contribution which the SP-CI 
attempts to compute very inefficiently. Hence clear ad- 
vantage in using URHF-CI versus SP-CI method. 

We now turn to the analysis of the excitation gaps. We 
focus on the energy gap between the spin singlet S — 
ground state and the spin triplet 5=1 excited state, the 
exchange energy J. 

In Figs. 1^1 the convergence of calculated exchange en- 
ergy J as a function of the number of configurations 
Mc/Nc is shown for increasing size Ng of single particle 
basis. We see that increasing the size Mc of the effective 
Hamiltonian initially leads to rapid oscillations in J fol- 
lowed by a smooth dependence. These calculations allow 
us to adjust Mc to extract numerically stable exchange 
energy J for each size of the single particle basis Ng. 

This allows us to compare the dependence of calcu- 
lated exchange energy using the URHF-CI and SP-CI 
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FIG. 11: Shown SP-CI and URHF-CI ground state ener- 
gies vs. number of single particle levels used in CI calcula- 
tion (Ns) for two electrons (top) and 4 electrons (bottom) at 
hu!c = lARy*. For two electrons, although, the ground state 
energy calculated from URHF is lower than single configu- 
ration SP-CI (top-left), but as the number of configurations 
increases, SP-CI method gives lower energy ground state (top- 
right). For four electrons, the URHF-CI ground state energy 
is lower than the one calculated by SP-CI for small number of 
configurations. As the number of configurations increases, the 
ground state energy calculated by SP-CI becomes the lowest. 



methods. Fig. ^| shows clearly fluctuations of J, cal- 
culated using SP-CI method , with increasing number of 
configurations. These fluctuations can be traced back to 
many level crossings in single particle orbitals. We are 
unable to extract reliable value of J using the commonly 
used SP-CI. By contrast, exchange energy J calculated 
using URHF-CI shows a smooth and convergent behavior 
as a function of Ng. 



HF-CI vs SP-CI - Dependence on the number 
of electrons N 



As we discussed in last subsection, computing the spec- 
trum of the CI Hamiltonian requires diagonalization of 
large matrices. The size of CI Hamiltonian matrices can 
be optimized significantly by a judicious choice of basis. 

In preceding subsection, we introduced URHF states 
as a suitable basis for CI method. In this subsec- 
tion, we present a systematic comparison between the 
ground state energy calculated by URHF-CI and SP-CI 
methods to remark that URHF-CI is a superior method 
to deal with a system with large number of electrons. 
We examine this in quantum dot molecules and in sin- 
gle dot parabolic confining potentials. In quantum dot 
molecules, we calculate the ground state energy of a 
system consisting of two to ten electrons. To compare 
URHF-CI and SP-CI, one must study the behavior of 
the spectrum as a function of Ng, and iV;(> Ng). With 
increasing number of configurations, this comparison is 
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FIG. 12: A comparison between SP-CI and URHF-CI ground 
state energies for six (top), eight (middle) and ten electrons 
(bottom) at TujJc ~ lARy* is shown. In contrast to two and 
four electron numbers, URHF-CI ground state energy is lower 
than SP-CI ground state energy for all configurations. 



Number of Electrons 

FIG. 13: Shown a zero magnetic field comparison between 
GS energies calculated by SP-CI and URHF-CI methods using 
Ns = 10, A*'; — 28 for a single dot parabolic potential {ujo = 
0.55Ry*) with electron numbers between A*' = 4, and 7. The 
calculated ground state energy of URHF-CI (SP-CI) is lower 
for a system with large (small) number of particles. 



carried out up to the point that ~ Ni, where all pos- 
sible HF configurations are being exhausted. In the limit 
of Ns = Ni, there exists a unitary transformation which 
maps URHF-CI Hamiltonian to SP-CI Hamiltonian, and 
thus the spectrum of SP-CI and URHF-CI become identi- 
cal. Because the number of configurations grows rapidly 
by increasing Ns, and because we would like to reach 
the limit of Ns — Ni, we construct URHF states out of 
small number of single particle levels. For the purpose 
of this comparison and without any loss of generality we 
present the results of our calculation using Ni = 8 in 
Figs. I11I12I In the case of two electrons, URHF energy 
is lower than the energy of SP-CI with single Slater deter- 
minant, Ns = 1, (see the top-left of Fig. Ill|) . By increas- 
ing the number of configurations we observe that SP- 
CI quickly lowers the ground state energy until Ns = 8 
where URHF-CI and SP-CI become equivalent. In the 
case of four electrons, the ground state energy of URHF- 
CI is lower for small number of configurations. Similar 
to the two electron system, with increasing number of 
configurations (beyond Ns = 5) SP-CI provides lower 
ground state energy. However, with increasing number 
of electrons from six up to ten, we find that URHF-CI 
method gives lower ground state energy within the whole 
range of Ns , as shown in Fig. 1121 

A similar comparison can be made for electrons in a 
single dot with parabolic confining potential. The energy 
difference calculated by SP-CI and URHF-CI, presented 
in Fig. ^1 reveals that for electron numbers between 
two and five, the ground state energy of SP-CI is lower 
than the ground state energy of URHF-CI. As expected, 
this order is reversed with increasing number of electrons 
(TV > 6). 

We therefore find URHF-CI method to work very well 



for quantum dot systems with large number of electrons. 

VI. CONCLUSION 

We have developed real space unrestricted hybrid 
Hartree Fock (HF) and configuration interaction (CI) 
method (URHF-CI) suitable for the calculation of ground 
and excited states of large number of electrons local- 
ized by complex gate potentials in quasi-two-dimensional 
quantum dot molecules. The effects of magnetic field and 
correlations are included in energies and directly in the 
many-particle wavefunctions making the method an at- 
tractive candidate for potential quantum information re- 
lated applications. The advantages of URHF-CI method 
over the commonly used CI method based on single par- 
ticle orbitals SP-CI are demonstrated. 
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VIII. APPENDIX: COULOMB MATRIX 
ELEMENTS 

In this appendix, we describe an efficient approach to 
calculate Coulomb matrix elements numerically. In the 
Hartree-Fock and configuration interaction method, the 
properties of the system are given by the single parti- 
cle spectrum and by the Coulomb matrix elements de- 
fined as two-electron integrals (see Eq. ^). For the 
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calculation involving 60 HF orbitals the total number of 
multi-dimensional integrals exceeds 12 x 10^. Here we 
describe an efBcicnt algorithm used in our calculation. 
In the numerical calculation of URHF-CI, the integrals 
of the Coulomb matrix elements are replaced by summa- 
tion over the grids r = {ihx,jhy,khz) (where {ijk} are 
integers, and {hx,hy,hz} are grid spacing) 

ri,r2 

^^T^ r,{r2)Mr2). (11) 

We further transform multi-dimensional wave function 
ipa{r) into a column vector Va(9) by mapping the 
multi-dimensional vector r onto a one- dimensional in- 



dex {ijk} q where q = 1,2,..., N^NyN^. Then the 
multi-dimensional integral can be converted into a vector- 
matrix multiplication 

Vc^p^, = J2 U{q, q') ^isM) (12) 

9.9' 

where ^av{q) = Va(9)V'i^(9) ^ vector containing all 
the possible pairs of single-particle wave functions, and 

2 

U is the matrix with elements times AriAfij. 

The sums in form of do-loops can be further parallelized. 
Due to the large dimension of matrix U, one can make 
use of domain decomposition to divide it into a number 
of smaller matrices and sum up the result of all the indi- 
vidual multiplications. 
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